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Abstract 

The  airborne  laser  (ABL)  uses  adaptive  optics  to  compensate  for  the  atmo¬ 
spheric  turbulence  between  the  ABL  and  the  target.  The  primary  purpose  of  this 
compensation  is  to  increase  the  energy  density  of  the  high  energy  laser  at  the  tar¬ 
get.  However,  the  specifics  of  the  engagement  scenario  require  the  tracking  point  of 
reference  and  the  adaptive  optics  point  of  reference  to  be  located  at  different  points 
on  the  target. 

This  research  considers  the  effects  of  tracking  a  target  in  one  direction  while 
compensating  for  atmospheric  turbulence  in  a  different  direction.  The  target  refer¬ 
ences  used  are  a  point  source  and  a  rectangle,  while  a  point  source  alone  is  used  for 
the  adaptive  optics  reference,  ft  will  be  shown  that  compensating  for  atmospheric 
turbulence  in  a  different  direction  than  tracking  results  in  a  bias  in  the  mean  tracking 
error  while  having  no  appreciable  affect  on  track  jitter. 


HIGH  ENERGY  LASER  POINTING 
THROUGH  EXTENDED  TURBULENCE 


I.  Introduction 

1.1  Problem,  Statement 

The  introduction  of  the  laser  brought  with  it  promises  of  Buck  Rogers  and 
Star  Wars.  Few  would  have  thought  that  the  largest  impact  would  have  been  in  how 
we  listen  to  music  or  watch  movies.  However,  developments  in  high  powered  lasers 
and  beam  control  have  revived  those  old  promises.  These  new  weapons,  such  as  the 
Airborne  Laser  (ABL),  seek  to  focus  megawatts  of  laser  energy  on  a  target  hundreds 
of  kilometers  away. 

Accomplishing  this  impressive  task  involves  solving  several  different  problems. 
The  first  and  most  obvious  entails  developing  a  laser  source  powerful  enough  to 
produce  megawatts  of  laser  energy.  After  that,  the  energy  needs  to  be  delivered 
accurately  to  the  target.  A  more  subtle  problem  involves  ensuring  that  the  laser 
energy,  when  it  reaches  the  target,  is  of  the  best  quality  possible  using  a  technique 
known  as  adaptive  optics.  Difficulties  may  arise  when  the  tracking  and  pointing 
becomes  dependent  on  the  performance  of  the  adaptive  optics.  In  the  case  of  the 
ABL,  the  tracking  reference  (i.e.,  the  nose  of  the  missile)  is  separated  from  the 
adaptive  optics  reference.  This  spatial  separation  causes  light  from  both  sources  to 
travel  different  paths.  The  purpose  of  this  research  is  to  determine  the  effect  on 
tracking  using  an  adaptive  optics  reference  that  is  not  co-located  with  the  tracking 
reference. 
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1.2  Scope 

Previous  work  [6,  7]  examined  the  effects  of  spatially  separated  paths  on  track¬ 
ing  alone  or  adaptive  optics  as  a  whole.  This  work  is  extended  by  tracking  solely 
along  one  path  and  compensating  for  higher  order  distortions  along  another  path. 
Performance  is  determined  by  sending  a  scoring  beam  back  through  the  atmosphere 
to  the  target.  The  adaptive  optics  algorithm  used  is  a  straightforward  least  squares 
reconstructor  [14],  which  provides  the  necessary  contrast  to  an  uncompensated  im¬ 
age. 

1.3  Thesis  Organization 

Chapter  two  contains  background  information  on  the  ABL,  adaptive  optics, 
and  image  tracking.  The  ABL  is  considered  first,  concentrating  on  the  geometry  of 
a  typical  ABL  engagement  and  the  optical  system  layout.  An  overview  of  adaptive 
optics  as  it  relates  to  this  research  follows.  Lastly,  a  few  tracking  techniques  are  bro¬ 
ken  down,  specifically  the  centroid  and  edge  tracking  algorithms.  Their  advantages 
and  disadvantages  are  discussed,  as  well  as  what  can  be  done  with  them. 

Chapter  three  looks  at  the  specific  methodology  used  for  determining  the  per¬ 
formance  of  the  systems.  First,  the  effects  of  anisoplanatic  compensation  are  exam¬ 
ined  theoretically.  Then,  the  simulation  itself  is  discussed  in  all  of  its  aspects. 

Chapter  four  examines  and  interprets  the  results  of  the  simulation  work.  This 
includes  validation  of  the  technique  as  well  as  results  from  the  tracking  experiments 
themselves.  It  also  includes  the  techniques  used  to  examine  the  results  and  how  to 
interpret  them. 

In  Chapter  five,  the  conclusions  are  discussed  and  avenues  for  future  research 
are  proposed. 
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II.  Background 


2. 1  Introduction 

Adaptive  optics  has  been  used  extensively  in  various  applications  including 
astronomical  imaging,  laser  beam  propagation,  and  even  within  the  laser  resonators 
themselves  [3,  5,  14,  15].  This  thesis  looks  at  a  particular  form  of  the  anisoplanatic  [6] 
problem  of  adaptive  optics  and  how  it  relates  to  the  Airborne  Laser  (ABL).  The 
discussion  begins  with  a  description  of  the  ABL  engagement  scenario  and  system 
layout.  This  is  followed  by  a  brief  overview  of  adaptive  optics.  This  will  include 
a  general  description  of  adaptive  optics  architecture,  some  limitations  imposed  by 
the  atmosphere,  and  how  tracking  is  separated  from  the  higher  order  correction.  A 
discussion  of  various  tracking  algorithms  follows.  Finally  a  description  of  the  specific 
problem  at  hand  will  tie  together  the  previous  sections  of  this  chapter. 

2.2  The  Airborne  Laser 

The  ABL  is  an  aircraft  designed  for  short  range  ballistic  missile  defense.  It 
consists  of  a  chemical  oxygen  iodine  laser  (COIL)  inside  a  Boeing  747  with  a  1.5 
meter  diameter  turret  mounted  on  the  nose  (Figure  2.1). 


Figure  2.1.  Artistic  rendition  of  the  Airborne  Laser. 


In  addition  to  the  COIL,  three  other  lasers  are  critical  to  system  operation.  The 
first  is  the  Active  Ranger  System,  which  provides  precise  range  data  for  tracking. 


2-1 


The  next  laser  is  the  track  illuminator  laser  (TILL).  The  TILL  flood  illuminates 
the  nose  of  the  missile  with  enough  energy  to  produce  a  return.  The  returning  light 
allows  the  ABL  to  determine  the  precise  position  of  the  nose  of  the  missile,  providing 
a  reference  from  which  all  of  the  remaining  lasers  are  aimed.  The  third  laser  is  the 
beacon  illuminator  laser  (BILL).  The  BILL  is  aimed  down  the  body  of  the  missile 
and  provides  the  reference  for  the  adaptive  optics.  Lastly  the  COIL  is  located  further 
down  the  body  of  the  missile.  The  entire  setup  is  shown  in  Figure  2.2. 


Figure  2.2.  ABL  engagement  scenario. 


The  separation  between  the  placement  of  the  BILL  on  the  missile  and  the 
placement  of  the  high  energy  laser  (HEL)  COIL  is  fixed  at  (2R/c  +  At)  v.  Here 
R  is  the  range  between  the  ABL  and  the  target,  c  is  the  speed  of  light,  At  is  the 
sum  total  of  any  latencies  in  the  ABL  beam  control  system,  and  v  is  the  transverse 
velocity  of  the  missile  relative  to  the  ABL.  This  distance  is  fixed  so  that  the  HEL 
travels  through  the  exact  same  path  that  the  reflection  of  the  BILL  from  the  missile 
travelled.  Therefore,  in  order  to  hit  a  desired  spot  on  the  missile,  the  BILL  aimpoint 
must  be  adjusted  [2], 

It  is  also  important  to  note  that  the  TILL,  BILL,  and  HEL  all  pass  through 
the  same  aperture  and  beam  control  system.  This  means  that  any  correction  applied 
as  a  result  of  the  tracking  system,  which  is  based  on  TILL  return,  effects  the  BILL 
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and  HEL.  Also,  any  adaptive  optics  correction  applied,  which  is  based  on  the  BILL 
return,  effects  the  TILL  and  HEL  also. 

2.3  Adaptive  Optics  Overview 

Adaptive  optics  involves  compensating  a  beam  of  light  for  distortions  that 
may  have  occurred  during  propagation.  Using  some  sort  of  reference,  ideally  a  point 
source  located  at  the  object  to  be  imaged,  the  distortions  experienced  by  the  reference 
can  be  compared  to  what  the  reference  would  have  looked  like  without  the  distortions. 
Once  these  differences  are  known,  action  can  be  taken  to  reverse  them.  Assuming 
light  from  the  object  travels  the  same  path  in  the  same  way  as  the  reference,  these 
actions  improve  the  quality  of  the  image  of  the  object. 

2.3.1  Adaptive  Optics  Architecture.  An  adaptive  optics  systems  can  be 
broken  into  three  very  basic  elements:  the  wavefront  sensor,  the  reconstructor,  and 
the  wavefront  corrector. 

1.  Wavefront  Sensor 

The  wavefront  sensor  is  the  device  used  to  measure  the  optical  phase  of  the 
light  received  from  the  reference.  There  are  many  ways  to  measure  the  incident 
phase  at  varying  degrees  of  complexity.  Wavefront  sensors  can  be  divided  into 
two  basic  classes:  image  plane  sensors  and  aperture  plane  sensors  [9]. 

Image  plane  sensors,  or  single  intensity  detectors,  attempt  to  maximize  the 
intensity  of  the  image  of  the  reference.  The  simplest  method  uses  a  single 
detector  to  measure  an  image  sharpness  function  such  as: 

S  =  Jl(x,y)M(x,y)dxdy  (2.1) 

where  I(x,y)  is  the  reference  intensity  at  the  image  plane  and  M(x,y)  is  some 
masking  function  such  as  a  pinhole.  The  detector  then  measures  the  total  light 


2-3 


passing  through  the  mask,  which  is  maximized  when  the  optical  system  is  free 
from  aberration  [9]. 

Aperture  plane  sensors  attempt  to  measure  the  wavefront  of  the  incoming 
beam.  There  are  several  methods  to  do  this,  including  Hartmann  sensors  and 
shearing  interferometers  [9].  A  Hartmann  wavefront  sensor  images  the  pupil 
of  the  optical  system.  By  using  a  lenslct  array  of  n  by  n  lenses,  the  pupil 
is  broken  up  into  n2  small  images,  or  subapertures,  on  a  detector  (shown  in 
Figure  2.3).  The  wavefront  sensor  then  calculates  the  position  of  the  spot 
in  each  subaperture,  and  uses  the  resulting  slopes  to  calculate  the  wavefront 
gradient  which  is  used  to  reconstruct  the  wavefront. 


Figure  2.3.  Hartmann  wavefront  sensor  showing  the  image  seen  at  the  detector 
plane,  an  array  of  spots  frorned  by  the  lenslet  array  [10]. 

2.  The  Reconstructor 

Once  the  measurements  are  made  by  the  wavefront  sensor,  the  actual  wavefront 
must  be  reconstructed.  This  is  done  by  the  wavefront  reconstructor.  Once  com¬ 
mon  way  to  reconstruct  the  wavefront,  and  the  method  used  in  this  research, 
is  to  use  a  least  squares  reconstructor  [14].  The  least  squares  reconstruction 
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matrix  is  defined  by 


M  =  ( HtH )  1  Ht 


(2.2) 


H  is  the  Jacobian  matrix  such  that  the  equation 

sm  H  c 


(2.3) 


is  satisfied.  Here  sm  is  a  column  vector  containing  the  measured  slopes  from 
the  wavefront  sensor,  and  c  is  a  column  vector  of  control  signals.  An  important 
note  on  the  least  square  reconstructor  is  that  neither  atmospheric  statistics  nor 
noise  statistics  are  required.  This  simplicity  makes  it  the  most  used  technique 
in  current  adaptive  optics  systems  [14]. 

3.  Wavefront  Corrector 

The  wavefront  corrector  is  a  device  that  can  alter  the  wavefront  characteristics 
of  the  incoming  wave.  Correctors  can  be  classed  as  either  modal  or  zonal. 
Modal  correctors  correct  specific  abberations.  These  abberations  are  typically 
described  using  the  Zernike  polynomials  [12]  defined  by 


Zevenj  =  y/n  +  1 R™  (r)  y/2  cos  rnO  I 

>  m  f  0 

Zoddj  =  y/n  +  1 R™  (r)  y/2  sin  md  J 

Zj  =  y/n  +  1  RQn  (r)  m  —  0 


(2.4) 


where 


-pm  _ 

rCn 


(n—m)/ 2 

E 

s=0 


(-l)s(n-2)! 


„n—2s 


s!  [(n  *f  m)/2  —  s]!  [(n  —  m)/2  —  s]! 


(2.5) 


The  values  of  n  and  m  are  always  integers  and  satisfy  m  <  n  and  n  —  \m\  = 
even.  The  index  j  is  used  to  denote  mode  ordering  and  is  a  function  of  n  and  m. 
The  Zernike  polynomials  are  typically  used  in  the  expansion  of  the  wavefront 
over  a  circular  aperture  of  radius  R.  If  (j)  (r,  6)  represents  the  wavefront  then 
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the  polynomial  expansion  over  a  circle  of  radius  R  is  given  by 


0  »)  =  aJ2J  (^’  <2J) 

3 

where  p  =  r/R  and  the  coefhcents  a3  are  given  by 

aj  —  f  d pW(p)<l>(Rp,e)Zj(p,e)  (2.7) 

where  W  (p)  is  the  pupil  function.  This  allows  the  wavefront  to  be  described  in 
terms  of  the  Zernike  polynomials,  for  instance  the  X  tilt  and  Y  tilt  components 
are  defined  by  the  Zernike  polymonials  Z 2  =  2 r  cos  9  and  Z3  =  2 r  sin  9  respec¬ 
tively,  defocus  described  by  ZA  =  (2r2  —  1),  and  so  on.  Zx  —  1,  referred  to 

as  piston,  is  unobservable  and  has  no  effect  on  the  image.  A  modal  corrector 
then  only  compensates  for  particular  abberations,  such  as  tilt  and  defocus. 

Zonal  correctors  break  the  wavefront  into  pieces  and  apply  the  necessary  cor¬ 
rection  to  each  piece  [15].  A  zonal  corrector  works  well  with  a  Hartmann 
wavefront  sensor,  in  that  each  subaperture  can  be  fitted  to  each  zone  of  the 
zonal  corrector. 

2.3.2  Atmospheric  Limitations.  Light  passing  through  any  inhomogeneous 
medium  becomes  distorted.  In  the  case  of  the  atmosphere,  these  distortions  are 
caused  by  minute  variations  in  the  index  of  refraction  along  the  beam’s  path.  These 
variations  are  ultimately  caused  by  temperature  fluctuations  resulting  in  density 
changes. 

The  turbulence  characteristics  of  a  given  path  through  the  atmosphere  is  com¬ 
monly  characterized  by  the  parameters  r0  and  90.  The  parameter  r0  is  known  as  the 
atmospheric  coherence  length  or  Fried  parameter  [6] .  This  number  gives  an  effective 
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diameter  of  the  system  if  there  were  no  atmosphere  present.  The  expression  for  r$  is 


r0  =  (2.905/6.88)  k 


dzC 2 


N 


»  [1  -  (z/L)f3 


-3/5 


(2,8) 


Here  k  is  the  wave  number  defined  by  2tt/X,  L  is  the  path  length,  is  the  index 
of  refraction  structure  constant,  and  z  is  the  integration  variable  along  the  path. 
Assumptions  are  that  there  is  spherical  wave  propagation  and  the  turbulence  follows 
the  Kolmogorov  spectrum  [15]  defined  by 


(. K )  =  0.033C/Ji~11/3 


(2,9) 


where  is  the  power  spectral  density  and  K  is  the  three-dimensional  spatial  wave 
number. 

The  parameter  9q  defines  the  isoplanatic  angle.  This  angle  determines  the 
angle  within  which  the  adaptive  optics  beacon  and  the  object  to  be  imaged  must 
exist  relative  to  the  observer,  outside  of  which  the  paths  are  uncorrelated.  Once  the 
beacon  and  object  are  separated  by  more  than  this  angle,  the  corrections  applied  by 
the  adaptive  optics  systems  begins  to  have  an  overall  degrading  effect  on  the  image 
of  the  object.  The  parameter  60  is  described  by  the  equation 

f  fL  }  -3/5 

60  =  l  2.905k2  J  dzC2N(z)z5/3i  (2.10) 

and  makes  the  same  assumptions  as  r o  [6]. 

2.3.3  Tracking  as  a  Separate  Problem.  Over  86%  of  the  total  wavefront 
error  is  due  to  the  tilt  component  of  the  atmosphere  [14].  Since  wavefront  correctors, 
such  as  deformable  mirrors,  have  small  dynamic  range,  tilt  is  usually  corrected  by 
using  a  fast  steering  mirror  (FSM).  Since  tilt  can  be  measured  from  a  single  image 
using  the  entire  aperture,  the  problem  boils  down  to  simply  tracking  the  center  of 
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the  reference,  assuming  the  reference  is  a  point  source.  In  most  cases,  the  object 
to  be  imaged  is  not  a  point  source,  but  some  other  adaptive  optics  reference  exists, 
or  can  be  made,  within  the  isoplanatic  angle  of  the  object.  In  these  instances,  it  is 
desirable  to  track  the  object  and  not  the  reference. 

When  tilt  is  corrected  using  a  separate  system,  it  becomes  necessary  not  to 
correct  for  tilt  using  the  adaptive  optics  portion  of  the  system  (i.e.,  the  wavefront 
sensor  and  corrector).  When  this  is  done,  the  wavefront  is  reconstructed  from  the 
wavefront  sensor  with  no  overall  X  or  Y  tilt.  This  is  also  know  as  tilt-removed 
reconstruction. 

2.4  Tracking  Algorithms 

Tracking  algorithms  can  be  broken  down  into  two  basic  forms,  centroid  trackers 
and  feature  trackers.  The  following  equations  are  the  definition  of  a  centroid: 


x  —  1  j  I(x,  y)  x  dx  dy  / 

!  j  J  I (x,y)  dx  dy 

(2.11) 

V=ff  I(x,y)ydxdy  y 

/  J  1 1(x,y)  dxdy 

(2.12) 

where  /  is  the  intensity  at  the  image  plane  of  the  receiver.  Feature  trackers  attempt 
to  match  the  image  to  some  reference  map.  This  technique  ranges  in  complexity 
from  a  simple  edge  tracker  to  a  full  blown  correlation  tracker. 

2.5  The  Problem 

The  specific  problem  examined  looks  at  the  implications  of  tracking  an  object 
in  one  direction,  while  compensating  for  everything  but  tilt  in  another  direction. 
Figure  2.4  shows  a  skeleton  of  the  simulation  experiment.  P\  and  P2  represent  point 
sources,  the  tracking  beacon  and  the  adaptive  optics  beacon  respectively,  at  the 
target  plane,  which  are  separated  by  Ax.  For  the  extended  source  simulations,  the 
tracking  beacon  is  replaced  by  an  extended  beacon  with  a  center  of  mass  located  at 
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Figure  2.4.  Schematic  of  the  basic  problem.  The  target  represents  two  point 
sources  at  different  locations  at  the  target  plane.  For  the  extended 
source  simulations,  one  of  the  point  sources  is  replaced  with  an  ex¬ 
tended  object. 

Pi  such  that  Ax  is  the  difference  between  the  center  of  the  target  and  the  adaptive 
optics  beacon. 

The  simulation  uses  a  Hartmann  wavefront  sensor  with  128  x  128  subapertures 
and  a  least  squares  estimator  with  the  tilt  components  removed.  The  wavefront  cor¬ 
rector  uses  an  optical  phase  delay  at  the  same  grid  level  as  the  wavefront  sensor  [16] . 
A  centroid  tracker  is  placed  at  the  image  plane  for  tilt  control.  Performance  is  de¬ 
termined  by  compensating  a  scoring  beam  using  the  tilt  commands  from  the  tracker 
and  sending  the  beam  to  a  target  board  (an  array  of  sensors  spaced  1cm  apart)  lo¬ 
cated  at  the  target  plane.  The  tracking  beacon,  whether  a  point  source  or  extended, 
is  always  located  at  the  center  of  the  target  board.  Further  details  will  be  discussed 
in  Chapter  3. 

2. 6  Summary 

This  chapter  introduced  several  principles  that  will  be  referenced  in  the  fol¬ 
lowing  chapters.  First,  the  operation  of  the  Airborne  Laser  was  discussed.  This  was 
followed  by  a  discussion  of  adaptive  optics.  The  adaptive  optics  section  reviewed 
the  basic  architecture  of  a  generic  system  and  some  basic  principles  of  operation. 
This  was  followed  by  a  discussion  on  the  limitations  that  the  atmosphere  imposes 
on  adaptive  optics  systems.  Next,  tracking  was  discussed,  focusing  on  the  centroid 
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tracker  and  how  it  is  used.  Finally,  a  discussion  of  the  problem  to  be  researched  was 
presented. 


2-10 


III.  Methodology 


3. 1  Introduction 

As  discussed  before,  this  research  investigates  a  particular  aspect  of  anisopla- 
natism.  The  experiment  is  done  entirely  using  wave  optics  simulation. 


3.2 

Analysis 

Beginning  with  the  equations  for  r0  and  0$,  these 

assuming  a  constant  altitude  resulting  in  a  constant  Cjj. 

can  be  further  simplified 

The  results  from  Fried  [6] 

are 

r0  =  3.02/c-6/5L-3/5  {C2n)~3/5 

(3.1) 

and 

d0  =  0.950Ar6/5L“8/5  (C2n)~3/5 

(3.2) 

k  is  the  wave  number  and  L  is  the  path  length.  The  simulation  is  actually  given  ro 
to  generate  the  turbulence  so  that 


0o  = 


0.950  r0 
3.02  T 


(3.3) 


The  simulation  uses  a  constant  altitude  for  two  primary  reasons.  First,  it  simplifies 
the  relationship  between  the  turbulence  level  in  the  simulation,  i.e.  r0,  and  the 
isoplanatic  angle.  Second,  it  closely  resembles  a  real  world  experimental  setup  at 
the  North  Oscura  Peak  (NOP)  facility  at  White  Sands  Missile  Range  [1].  This 
facility  consists  of  a  1.0  meter  (m)  telescope,  with  associated  tracking  and  adaptive 
optics  systems,  located  on  a  mountain  peak  at  an  altitude  of  approximately  2400m.  A 
second  station  at  Salinas  peak,  nearly  50  kilometers  (km)  away  provides  a  stationary 
target,  which  can  be  used  by  the  facilities  at  NOP  for  tests. 
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3. 3  Simulation  Method 


The  simulation  used  the  wave  optics  package  Wave  Train™  [11],  The  simula¬ 
tion  uses  two  dimensional  meshes  of  complex  numbers  to  represent  the  propagating 
bearn(s)  of  light  [4],  For  propagation  of  the  complex  meshes,  WaveTrain  uses  the 
Fast  Fourier  Transform  technique  base  on  the  Fresnel  diffraction  integral  [8]: 

OO  OO 

U  (x,  y )  =  J  J  {u  (f ,  V)  ej^eW) }  e~j^+y^d^dri  (3.4) 

— OO  — OO 

where:  x  and  y  represent  the  coordinates  at  the  output  plane,  C,  and  rj  represent  the 
coordinates  at  the  input  plane,  and  z  is  the  propagation  distance. 

The  simulation  used  is  comprised  of  three  major  subsystems:  the  target,  the 
atmosphere,  and  the  observer,  as  shown  in  Figure  3.1.  Each  subsystem  contains 
various  components  of  the  simulation  such  as  sources,  detectors,  and  other  optical 
subsystems.  Figure  3.1  also  shows  the  lines  connecting  the  subsystems  together, 
representing  which  components  share  information  and  in  what  direction.  In  this 
case,  each  connection  shows  a  two-way  relationship,  meaning  data  is  sent  in  both 
directions  along  the  connection.  Also,  blocks  are  given  two  names.  The  first  name, 
which  is  always  capitalized,  describes  the  block’s  C++  class.  The  second  name,  in  all 
lowercase  letters,  defines  an  instance  of  the  class  that  defines  the  block.  Wavetrain 
uses  the  graphical  interface  to  define  relationships  between  subsystems,  which  is 
used  to  then  generate  C++  source  code,  which  is  compiled  and  executed  for  each 
simulation  run. 

3.3.1  The  Target.  The  target  subsystem,  shown  in  Figure  3.2,  consists  of 
the  target  board,  the  target,  and  the  adaptive  optics  reference.  The  large  arrows 
in  Figure  3.2  represent  the  light  coming  into  the  target  subsystem.  The  labels  in¬ 
coming  Transmitted  and  outgoinglncident  represent  the  light  leaving  the  target  block 
and  entering  the  target  block  respectively.  Incoming  is  defined  as  light  moving  from 
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Figure  3.1.  Overall  system  layout. 

right  to  left,  and  outgoing  is  defined  as  light  moving  from  left  to  right.  An  easy  way 
to  remember  is  that  everything  is  defined  from  the  point  of  view  of  the  HEL  on  the 
ABL. 

Starting  with  the  outgoing  Incident  arrow,  the  light  is  sent  to  an  Incoming- 
Combiner  block.  This  block  allows  for  two  sources  to  be  combined  for  transmittance 
in  one  direction,  while  the  other  direction  just  allows  one  beam  to  pass.  The  light 
from  the  outgoing  Incident  arrow  passes  through  the  Incoming  Combiner  block  un¬ 
modified  and  is  sent  to  the  target  board.  From  the  other  direction,  the  light  from 
the  PointSource  block  and  the  incoherentreflector  block  are  combined  and  sent  to 
the  incoming  Transmitted  arrow  and  out  the  Target  subsystem. 

3. 3. 1.1  Target  Board.  The  target  board  consists  of  a  grid  of  256  by 
256  sensors,  which  will  be  the  propagation  grid  of  the  simulation.  Each  sensor  is  1  cm 
wide  and  is  sensitive  only  light  to  with  a  wavelength  of  1  /jrn,  the  wavelength  of  the 
HEL.  1  /irn  was  used  as  the  center  wavelength  of  the  simulation  since  it  is  the  same 
wavelength  used  at  the  NOP  facility.  The  target  board  feeds  a  centroid  calculator, 
the  results  of  which  are  stored  for  post-processing.  Additionally,  the  images  of  the 
target  board  can  be  recorded. 

3.3. 1.2  The  Target.  The  target  shown  in  Figure  3.2  consists  of  an 
incoherent  reflector  and  a  uniform  illuminator.  The  reflector  consists  of  a  32  pixel  tall 


TargetBoard 


Centroid 


Figure  3.2.  Target  subsystem  layout. 

by  16  pixel  wide  array.  The  incoherent  reflector  acts  as  a  light  scattering  device,  with 
no  coherence  to  the  spatial  properties  of  light  reflected.  Eight  random  realizations 
were  used  to  produce  each  incoherent  reflection  pattern.  The  uniform  illuminator 
simply  flood  illuminates  the  incoherent  reflector  with  1.001  /irn  wavelength  light.  The 
array  is  centered  on  the  the  center  of  the  target  board.  For  the  point  source  tracking 
portion  of  the  experiment,  the  incoherent  reflector  and  the  uniform  illuminator  are 
replaced  by  a  point  source  centered  on  the  target  board.  This  takes  the  place  of  the 
TILL  on  the  ABL. 

3.3. 1.3  The  Adaptive  Optics  Reference.  This  experiment  used  a  point 
source  for  the  adaptive  optics  reference.  The  source  has  a  wavelength  of  0.999  /irn 
and  can  be  placed  anywhere  on  the  target  board.  This  experiment  holds  the  source 
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at  the  vertical  center  of  the  target  board  and  moves  along  the  horizontal  axis.  This 
takes  the  place  of  the  BILL  on  the  ABL. 

The  main  reason  different  wavelengths  were  chosen  for  the  sources  was  to  allow 
the  sensors  in  the  simulation  to  discriminate  between  them,  thus  only  a  0.001  fin i 
difference.  Also,  by  placing  all  of  the  sources  at  the  target  end  of  the  simulation, 
they  represent  ideal  illumination.  In  the  case  of  the  ABL,  the  BILL  and  TILL  are 
distorted  by  the  atmosphere  before  they  reach  the  target.  In  this  simulation  this  is 
not  the  case.  This  was  done  to  isolate  the  anisoplanatic  effects  from  any  that  might 
arise  due  to  the  outgoing  distortion  of  the  BILL  or  TILL. 

3.3.2  The  Atmosphere.  The  atmosphere  is  modelled  using  several  discrete 
thin  phase  screens  place  along  the  path  between  the  target  and  the  observer.  A  thin 
phase  screen  simulates  the  atmosphere  by  generating  a  grid  such  that 

U  (x,  y)  =  U0  (x,  y)  Ts  (x,  y)  (3.5) 

where  Uq  is  the  input  field,  U  is  the  output  field,  and  Ts  is  the  transparency  function 
describing  the  random  field  perturbation  due  to  the  atmosphere  [14].  Thin  phase 
screens  are  useful  in  wave  optics  simulations,  in  that  an  entire  random  segment  of 
atmosphere  can  be  reduced  to  one  simple  relation,  vastly  reducing  computation  time. 
This  particular  simulation  used  a  propagation  length  of  50km  with  5  phase  screens 
evenly  spaced  along  the  path.  The  number  of  phase  screens  was  chosen  based  on  the 
rule  of  thumb  that  dictates  one  phase  screen  for  every  10km  of  propation. 

Figure  3.3  shows  the  internals  of  the  Atmosphere  subsystem.  The  incoming 
light  follows  a  path  from  the  top  right  arrow,  labeled  WaveTrain,  to  the  bottom 
left  arrow.  The  outgoing  light  goes  from  the  top  left  arrow  to  the  bottom  right  ar¬ 
row.  The  TransverseVelocity  blocks  allow  for  the  addition  of  the  simulation  of  wind. 
It  does  this  by  “sliding”  the  phase  screens,  which  are  contained  in  the  AtmoPath 


3-5 


block,  across  the  path.  This  simulation  does  not  contain  any  wind  velocity,  thus  the 
atmosphere  is  “frozen”  for  the  duration  of  each  run. 


Figure  3.3.  Atmospheric  turbulence  block. 


3.3.3  The  Observer.  The  observer  block  consists  of  several  subcomponents. 
Overall,  it  consists  of  four  cameras  with  centroid  detectors,  the  adaptive  optics  com¬ 
pensation,  the  outgoing  steering  mirror,  and  the  scoring  laser.  Starting  from  the 
bottom  right  of  Figure  3.4,  the  light  enters  the  telescope  block.  The  telescope  has 
diameter  of  .5m  with  no  annulus  and  focused  on  the  target,  50km  away.  The  tele¬ 
scope  therefore  collimates  the  light  from  that  point  on  back.  The  incoming  light  then 
passes  through  the  Outgoing  Combiner  block  and  proceeds  to  a  beamsplitter.  This 
splits  the  beam  into  two  nearly  identical  paths,  the  difference  being  that  one  path  has 
a  compensator  and  the  other  does  not.  The  output  of  both  of  the  camera  blocks  are 
tilt  commands  that  are  sent  to  a  selector.  The  selector  determines  which  one  of  the 
four  cameras  will  control  the  beam  steering  mirror.  The  four  cameras  represent  the 
compensated  TILL  image,  uncompensated  TILL  image,  compensated  BILL  image, 
and  uncompensated  BILL  image.  Only  one  camera  is  allowed  to  control  the  steering 
mirror  at  a  time  so  that  comparisons  can  be  made  between  the  performance  of  the 
different  cameras.  The  selector  is  designed  so  that  tilt  commands  can  be  combined 
between  the  cameras  if  desired.  The  scoring  laser  (shown  as  the  bottom  left  block) 
is  sent  to  the  steering  mirror,  through  the  telescope,  and  out  to  the  target  board. 

One  important  aspect  of  the  observer  is  the  limitation  of  the  steering  mirror. 
The  simulation  allows  for  the  steering  mirror  to  have  an  instantaneous  response  to  the 
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commanded  input.  However,  the  light  takes  one  time  step  to  propagate.  Therefore, 
in  the  first  time  step,  the  light  propagates  from  the  target  to  the  observer.  In  the 
next  time  step,  the  cameras  calculate  the  centroids  and  the  command  is  sent  to  the 
steering  mirror.  In  the  third  time  step  the  mirror  reaches  its  final  position,  taking 
only  the  one  time  step  to  do  so.  Finally,  in  the  fourth  time  step  the  light  from  the 
HEL  is  propagated  to  the  target  board,  after  having  been  redirected  by  the  steering 
mirror.  In  all,  4  time  steps  are  required,  not  for  the  simulation  of  wind  or  any  other 
temporal  effects  of  the  atmosphere,  but  to  allow  the  control  elements  to  reach  their 
final  state  and  for  the  light  to  propagate. 

3.3.3. 1  The  Compensator.  Figure  3.5  shows  the  adaptive  optics 
compensation  block.  Following  the  light  in  from  the  right,  the  incoming  beam  is 
split.  The  first  path,  to  the  left,  is  again  split  between  a  beam  steering  mirror 
and  a  camera.  The  camera  is  sensitive  only  to  the  wavelength  of  the  the  adaptive 
optics  beacon.  The  camera  outputs  an  intensity  map  to  the  centroid  calculator, 
which  outputs  x  and  y  tilts.  These  tilts  are  inverted  and  sent  as  commands  to 
the  BeamSteering Mirror  block,  with  a  time  constant  of  1  /is.  This  process  removes 
the  tilt  from  the  beam  coming  from  the  second  beam  splitter,  thus  producing  a 
tilt-removed  reconstructor  while  still  using  a  simple  least  squares  reconstructor. 

After  the  tilt  has  been  removed,  the  field  is  sent  to  a  SimpleFieldSensor  block. 
This  block,  just  as  the  name  implies,  senses  the  field  at  the  grid  level.  The  resulting 
values  are  sent  to  the  reconstructor  block.  This  reconstructor  was  written  entirely  in 
Matlab®  by  Jeff  Barchers  and  Brent  Ellcrbrook  from  the  Air  Force  Research  Labo¬ 
ratory,  Starfire  Optical  Range  and  modified  for  WaveTrain  by  Virgil  Zetterlind  [16]. 
It  implements  the  least  squares  algorithm  discussed  in  Chapter  2.  The  reconstructor 
outputs  commands  to  an  optical  phase  delay  map,  that  acts  on  the  original  incoming 
beams  (both  the  TILL  and  BILL  returns),  thereby  compensating  the  beam. 
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Figure  3.5.  Compensator  block. 

3. 3. 3. 2  The  Cameras.  The  camera  blocks,  shown  in  Figure  3.6  and 
labeled  comp  and  uncomp  in  Figure  3.4  ,  are  identical  to  each  other  with  the  excep¬ 
tion  of  their  inputs.  One  receives  the  beam  unmodified  after  the  telescope  while  the 
other  receives  the  beam  after  it  has  been  compensated  for  higher  order  (other  than 
tilt)  aberrations.  In  both  camera  blocks,  the  top  sensor  is  tuned  only  to  the  track 
beacon,  while  the  bottom  is  tuned  to  the  adaptive  optics  beacon.  The  outputs  of 
the  subsystem  are  the  centroids  from  the  two  sensors. 

3. 3. 3. 3  The  Control  Selector.  The  last  subsystem  within  the  observer 
is  the  control  selector  block.  In  this  block,  the  centroid  calculation  from  the  two 
camera  blocks  (totalling  four  in  all:  uncompensated  track  beacon,  uncompensated 
adaptive  optics  beacon,  compensated  track  beacon,  and  compensated  adaptive  optics 
beacon)  are  multiplied  by  some  gain  and  then  summed  together.  By  toggling  the  gain 
at  runtime  between  0  and  1,  the  camera  that  controls  the  scoring  laser  steering  mirror 
can  be  chosen.  For  example,  by  setting  the  gains  to  {1,0, 0,0}  the  uncompensated 
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Figure  3.6.  Camera  blocks. 

track  camera  will  control  the  mirror,  and  by  using  {0,0, 1,0}  the  compensated  track 
camera  is  selected.  The  sign  inversion  necessary  for  proper  command  of  the  mirror 
is  inserted  in  the  last  summation  block  just  before  the  output  of  the  subsystem. 
This  technique  not  only  allows  for  the  selection  of  one  set  of  commands,  but  also 
allows  for  the  combination  of  commands  from  multiple  sources.  For  example,  by 
using  {0.5, 0,0. 5,0}  the  average  between  the  uncompensated  track  camera  and  the 
compensated  track  camera  can  be  used  to  command  the  steering  mirror. 

3-4  Summary 

In  this  chapter,  we  have  developed  some  relations  to  describe  the  turbulence 
levels  in  the  experiment.  We  have  also  described  the  particular  simulation  technique 
that  was  utilized.  We  reviewed  each  segment  of  the  simulation,  how  it  works,  and 
how  the  segments  were  set  up  to  work  together.  We  also  noted  that  each  random 
realization  was  4  time  steps  long,  allowing  the  delays  from  all  of  the  sensors  to  the 
control  elements  to  propagate  through  and  bring  the  system  to  a  steady  state.  The 
next  chapter  will  look  at  the  results  from  all  of  the  simulations. 
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Figure  3.7.  The  control  selector. 
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IV.  Results 


4-1  Introduction 

This  chapter  examines  the  results  of  the  propagation  simulations  that  were 
conducted.  It  is  divided  into  two  main  sections  dealing  with  the  point  source  ex¬ 
periments  first  and  then  the  extended  source  experiments.  The  point  source  and 
extended  source  were  done  in  the  same  manner  using  the  same  random  seeds.  Each 
set  consisted  of  80  random  samples,  21  different  separations  between  the  track  bea¬ 
con  and  adaptive  optics  beacon,  and  10  values  of  r0  for  compensated  tracking.  The 
complete  list  of  system  parameters  are  shown  in  Table  4.1.  The  uncompensated  sets 
were  identical  save  for  the  fact  that  only  one  separation  was  used,  since  the  place¬ 
ment  of  the  adaptive  optics  beacon  was  irrelevant  to  the  uncompensated  results. 


Parameter 

Value  (s) 

Units 

r0 

{10,9,8,7,6,5,4,3,2,1} 

cm 

No.  of  Realizations 

80 

NA 

Beacon  Separations 

{0,2.5, 5, 7.5,.. .,47.5,50} 

cm 

Range 

50 

km 

No.  of  Phase  Screens 

5 

NA 

HEL  Wavelength 

1.000 

/jrn 

BILL  Wavelength 

0.999 

/jrn 

TILL  Wavelength 

1.001 

/JH1 

FSM  Time  Constants 

1 

/is 

Target  Board  Resolution 

1 

cm 

Target  Board  Size  (H  x  W) 

2.56  x  2.56 

m 

Propagation  Grid  Size 

256  x  256 

NA 

Propagation  Grid  Spacing 

1 

cm 

Camera  IFOV 

0.2 

yurad 

Camera  Size 

256  x  256 

pixels 

Platform  Height 

2400 

m 

Target  Height 

2400 

m 

Table  4.1.  Parameters  used  for  the  simulation. 
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4-2  Turbulence  Levels 

As  it  turns  out,  the  turbulence  levels  chosen  for  these  simulation  runs  tend  to 
be  very  strong.  The  method  used  for  determining  the  turbulence  level  in  WaveTrain 
under-represents  the  actual  level  of  turbulence.  The  parameter  ro,  developed  origi¬ 
nally  for  astronomical  type  viewing,  takes  into  consideration  primarily  the  near  field 
turbulence.  For  a  vertical  path  this  is  fine,  but  for  horizontal  paths  it  is  not.  Table 
4.2  shows  the  chosen  values  of  ro  compared  with  other  methods  of  measuring  tur¬ 
bulence.  In  this  table,  along  side  r0  and  90,  the  corresponding  values  for  C„,  Rytov 


r o  (cm) 

90  (/irad) 

r‘2 

n 

Rytov  Variance 

Clear  1  Factor 

10 

0.63 

1.48xl0~16 

0.642 

2.97 

9 

0.57 

1.77xl0“16 

0.765 

3.54 

8 

0.50 

2.15AL0-16 

0.931 

4.30 

7 

0.44 

2.69a;10~16 

1.16 

5.38 

6 

0.38 

3.48xl0”16 

1.50 

6.95 

5 

0.31 

4.71xl0”16 

2.04 

9.42 

4 

0.25 

6.83a;10"16 

2.95 

13.7 

3 

0.19 

l.lOaTO”15 

4.77 

22.1 

2 

0.13 

2.17AL0-15 

9.38 

43.4 

1 

0.063 

6.89a;10”15 

29.8 

138 

Table  4.2.  Other  parameters  for  measuring  turbulence  level.  Even  at  the  weakest 
turbulence  level  used,  the  level  of  turbulence  is  very  high. 


variance,  and  Clear  1  Factor  are  shown.  was  derived  using  Equation  (  3.1).  The 
Rytov  variance  [3]  represents  the  irradiance  fluctuations  of  the  wave  at  the  receiver. 
It  is  defined  by  the  equation 


of  =  O.mC^L11/6 


(4.1) 


Weak  fluctuations  are  associated  with  af  <C  .1  while  strong  fluctuation  are  associated 
with  af  .1.  Table  4.2  indicates  that  the  turbulent  levels  chosen  are  clearly  in  the 
strong  fluctuation  region.  The  other  parameter,  Clear  1  Factor,  is  simply  a  scaling 
of  by  a  value  of  5xl0-1'.  A  Clear  1  Factor  of  1  is  the  baseline  design  criterion  for 
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the  Airborne  Laser.  Table  4.2  shows  that  the  turbulence  levels  chosen  go  far  beyond 
the  reasonable  conditions  that  the  Airborne  Laser  would  experience.  It  would  have 
been  better  to  determine  the  turbulence  strength  by  C %  and  Rytov  variance,  as  well 
as  r0. 


4-3  Point  Source  Results 

The  point  source  results  follow.  First  some  sample  images  are  shown,  and  then 
the  pointing  performance  data. 

4-3.1  Test  Images.  Figures  4.1  through  4.3  show  the  uncompensated  track 
beacon  images  at  the  observer.  Each  sub-figure  is  composed  of  four  random  samples 
at  the  same  turbulence  level.  Across  the  sub-figures  the  same  random  seed  is  shown. 
This  makes  it  appear  like  they  are  all  pictures  of  the  same  object,  just  progressively 
blurrier.  As  can  be  seen  in  Figure  4.3,  the  cases  in  which  r0  is  2cm  or  less,  the  beacon 
begins  to  have  portions  outside  the  viewing  area.  For  this  reason  the  data  at  those 
extreme  turbulence  levels  are  omitted  from  the  rest  of  the  results.  Also,  even  for 
the  weaker  turbulence  levels,  significant  spot  breakup  can  be  noted.  This  speckling 
of  the  image  causes  bright  and  dark  spots  to  appear,  which  could  have  a  significant 
effect  on  the  centroid,  since  it  is  an  intensity  weighted  measurement.  This  speckling 
would  also  affect  the  performance  of  the  adaptive  optics  compensation.  Figures  4.4 
and  4.5  show  the  compensated  images  of  the  track  beacon. 

Figures  4.6  and  4.7  show  the  scoring  beam  as  seen  by  the  target  board.  These 
are  shown  in  the  same  way  as  in  Figures  4.1  through  4.3,  with  four  random  samples 
per  sub-figure.  The  overall  impression  is  the  same,  with  the  ultimate  result  being 
that  data  at  r0  =  2 cm  and  r0  =  lcm  is  not  able  to  be  captured  accurately. 

4-3.2  Tilt  Performance.  Tracking  performance  is  often  measured  in  terms 
of  jitter,  or  tilt  variance.  Figures  4.8  and  4.9  show  resulting  error  variance  at  the 
target  board  for  both  compensated  and  uncompensated  tracking.  The  graphs  show 
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Figure  4.1.  Sample  images  of  the  point  source  track  beacon.  Each  sub- figure  shows 
four  different  realizations,  (a)  ro  =  10cm  (b)  r$  =  9cm  (c)  tq  =  8cm 
(d)  r0  =  7cm 
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Figure  4.2.  More  sample  images  of  the  point  source  track  beacon,  (a)  ro  =  6cm 
(b)  r0  =  5cm  (c)  r0  =  4cm  (d)  r0  =  3cm 
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Figure  4.3.  Sample  images  of  the  point 
lcm 


variance  as  a  function  of  separation  between  the  adaptive  optics  beacon  and  the 
target  beacon.  As  expected,  the  tracking  performance  increases  (tilt  variance  is 
reduced)  when  the  higher  order  compensator  is  used.  There  are  a  couple  of  notable 
abnormalities  in  the  results  for  tilt  variance.  These  occur  on  Figure  4.8  (a)  and 
(c).  These  spikes,  in  each  of  the  three  cases,  result  from  a  single  run  producing  an 
abnormally  high  tilt  error  (the  difference  between  the  center  of  the  target  board  and 
the  measured  centroid).  It  appears  that  in  each  case  the  steering  mirror  had  not 
settled,  causing  the  large  error.  The  steering  mirrors  all  have  time  constants  of  1/xs, 
meaning  each  would  reach  a  steady  state  value  before  the  next  time  step.  However, 
in  each  of  these  anomalies  the  tilt  error  is  exactly  the  same  for  the  third  and  fourth 
time  step  (remembering  that  each  run  is  four  time  steps  long).  This  would  indicate 
that  the  steering  mirror  did  not  receive  the  command  to  move  before  the  fourth  time 
step.  All  attempts  to  recreate  these  anomalies  failed. 

Figures  4.10  through  4.13  show  the  tilt  error  for  first  the  x  axis  and  then  the 
y  axis.  Keep  in  mind  the  adaptive  optics  beacon  was  moved  along  the  x  axis.  The 
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Figure  4.4.  Sample  images  of  the  compensated  point  source  track  beacon.  Note 
that  the  adaptive  optics  system  degradation  as  the  turbulence  level 
increased  (a)  r0  =  10cm  (b)  r0  =  9cm  (c)  r o  =  8cm  (d)  r0  =  7cm 
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Figure  4.5. 


More  sample  images  of  the  compensated  point  source  track  beacon, 
(a)  r0  =  6cm  (b)  r0  =  5  cm  (c)  r0  =  4cm  (d)  r0  =  3cm 
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Figure  4.6. 


(c) 

Sample  images  at  the  target  hoard,  (a)  r0 
r0  =  8cm  (d)  r0  =  7cm 


(d) 

10cm  (h)  r  o  =  9cm  (c) 
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(c)  (d) 

Figure  4.7.  Sample  images  at  the  target  board,  (a)  r0  =  6cm  (b)  r0  =  5cm  (c) 
r0  =  4cm  (d)  r0  =  3cm 
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(c) 
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Figure  4.8.  Tilt  variance.  Each  sub-figure  shows  the  tilt  variance  for  compensated 
and  uncompensated  (a  horizontal  line  since  it  is  independent  of  beacon 
separation )  as  a  function  of  separation  between  the  adaptive  optics 
beacon  and  the  track  beacon.  Beacon  separation  is  measured  in  values 
of  do-  The  spikes  in  (a)  and  (c)  are  the  results  of  a  single  run  in  each 
case.  These  anomalies  were  likely  the  results  of  the  steering  mirror  not 
reaching  its  final  state  properly.  Attempts  to  reproduce  the  specihc 
error  failed,  (a)  ro  =  10cm  (b)  r o  =  9cm  (c)  ro  =  8cm  (d)  Tq  =  7cm 
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Variance  (|irad  )  Variance  (|irad  ) 


Tilt  Variance  at  Target,  rQ  =  0.06 


Tilt  Variance  at  Target,  rQ  =  0.05 


Figure  4.9.  Tilt  variance.  Each  sub-figure  shows  the  tilt  variance  for  compensated 
and  uncompensated  (a  horizontal  line  since  it  is  independent  of  beacon 
separation )  as  a  function  of  separation  between  the  adaptive  optics 
beacon  and  the  track  beacon.  Beacon  separation  is  measured  in  values 
of  60.  (a)  r0  =  6cm  (b)  r0  =  5cm  (c)  r0  =  4cm  (d)  r0  =  3cm 
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dominant  feature  in  the  x  axis  is  the  tendency  to  move  towards  the  positive  as  the 
beacon  moves  more  negative.  The  likely  explanation  is  an  overcorrection  of  tilt  by 
the  tilt-removal  section  of  the  compensator  block.  This  is  probably  a  result  of  part 
of  the  return  from  the  beacon  being  moved  off  the  detector.  This  would  then  result 
in  an  under-representation  of  the  actual  tilt,  therefore  causing  the  error.  The  other 
notable  effect  is  the  apparent  bias  in  the  y  axis  that  appears  as  the  beacon  is  moved. 
It  seems  that  the  effect  is  a  result  of  the  anisoplanatic  compensation  and  not  from  an 
artifact  of  the  simulation.  However,  if  the  anisoplanatic  compensation  were  having 
an  effect  on  the  ability  of  the  system  to  track  the  target,  then  the  variance  would 
have  increased  as  well. 

4-4  Extended  Source  Results 

The  extended  source  results  follow.  First  some  track  beacon  and  target  board 
image  samples  are  shown,  and  then  the  pointing  performance  data. 

4-4-1  Test  Images.  Figures  4.14  and  4.15  show  the  uncompensated  track 
beacon  images  at  the  observer.  Each  sub-figure  is  actually  composed  of  four  random 
samples  at  the  same  turbulence  level.  Across  the  sub-figures  the  same  random  seed 
is  shown.  This  makes  it  appear  like  they  are  all  pictures  of  the  same  object,  just 
progressively  blurrier.  It  is  also  interesting  to  observe  the  track  beacon  using  com¬ 
pensated  images,  shown  in  Figures  4.16  and  4.17.  The  shape  of  the  32cm  by  16cm 
rectangle  target  can  be  clearly  seen  in  the  lower  turbulence  cases. 

4-4-2  Tilt  Performance.  Figures  4.18  and  4.19  show  resulting  error  variance 
at  the  target  board  for  both  compensated  and  uncompensated  tracking.  Figures  4.20 
through  4.23  show  the  tilt  error  for  first  the  x  axis  and  then  the  y  axis. 

The  results  for  the  extended  source  are  very  similar  to  that  of  the  point  source. 
This  is  likely  due  to  the  method  of  tracking  and  the  nature  of  the  target.  In  both  cases 
a  centroid  tracker  was  used  to  track  a  symmetrical  target.  Without  the  tracker  being 
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(c) 


Compensated  X  Tilt  Error  at  Target  rQ  =  0.09 


(b) 

Compensated  X  Tilt  Error  at  Target  rQ  =  0.07 


(d) 


Figure  4.10.  X  Tilt  error.  Each  sub-figure  shows  the  tilt  error  for  compensated  and 
uncompensated  (a  horizontal  line  since  it  is  independent  of  beacon 
separation )  as  a  function  of  separation  between  the  adaptive  optics 
beacon  and  the  track  beacon.  Beacon  separation  is  measured  in  values 
of  9q.  (a)  r  o  =  10cm  (b)  r  o  =  9cm  (c)  r  o  =  8cm  (d)  r0  =  7cm 
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Error  (m)  Error  (m) 


Compensated  X  Tilt  Error  at  Target,  r0  =  0.06 


(a) 

Compensated  X  Tilt  Error  at  Target,  r0  =  0.04 


(C) 


Compensated  X  Tilt  Error  at  Target  rQ  =  0.05 


(b) 

Compensated  X  Tilt  Error  at  Target  rQ  =  0.03 


(d) 


Figure  4.11.  Tilt  error.  Each  sub-figure  shows  the  tilt  error  for  compensated  and 
uncompensated  (a  horizontal  line  since  it  is  independent  of  beacon 
separation )  as  a  function  of  separation  between  the  adaptive  optics 
beacon  and  the  track  beacon.  Beacon  separation  is  measured  in  values 
of  0O-  (&)  f o  =  6cm  (b)  r0  =  5cm  (c)  r0  =  4cm  (d)  r0  =  3cm 
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Compensated  Y  Tilt  Error  at  Target  rQ  =  0.1  Compensated  Y  Tilt  Error  at  Target  r0  =  0.09 
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Figure  4.12.  Y  Tilt  error.  Each  sub-figure  shows  the  tilt  error  for  compensated  and 
uncompensated  (a  horizontal  line  since  it  is  independent  of  beacon 
separation )  as  a  function  of  separation  between  the  adaptive  optics 
beacon  and  the  track  beacon.  Beacon  separation  is  measured  in  values 
of  d0  ■  (&)  fo  =  10cm  (b)  r  o  =  9cm  (c)  r0  =  8cm  (d)  r0  =  7cm 


4-16 


Error  (m)  Error  (m) 


Compensated  Y  Tilt  Error  at  Target,  r0  =  0.06 
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Compensated  Y  Tilt  Error  at  Target,  r0  =  0.04 
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(C) 


Compensated  Y  Tilt  Error  at  Target  rQ  =  0.05 


0  5  10  15 

Separation  (0O) 


(b) 


Compensated  Y  Tilt  Error  at  Target  rQ  =  0.03 
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(d) 


Figure  4.13.  Y  Tilt  error.  Each  sub-figure  shows  the  tilt  error  for  compensated  and 
uncompensated  (a  horizontal  line  since  it  is  independent  of  beacon 
separation )  as  a  function  of  separation  between  the  adaptive  optics 
beacon  and  the  track  beacon.  Beacon  separation  is  measured  in  values 
of  0O-  (&)  f o  =  6cm  (b)  r0  =  5cm  (c)  r0  =  4cm  (d)  r0  =  3cm 
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(c)  (d) 

Figure  4.14.  Sample  images  of  the  uncompensated  extended  source  track  beacon. 

Each  sub- figure  shows  four  different  realizations,  (a)  r$  =  10cm  (b) 
r0  =  9cm  (c)  r0  =  8cm  (d)  r0  =  7cm 
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Figure  4.15. 


Sample  images  of  the  uncompensated  extended  source  track  beacon, 
(a)  r0  =  6cm  (b)  r0  =  5cm  (c)  r0  =  4cm  (d)  r0  =  3cm 
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Figure  4.17.  Sample  images  of  the  compensated  extended  source  track  beacon,  (a) 
r0  =  6cm  (b)  Tq  =  5  cm  (c)  r0  =  4cm  (d)  r0  =  3cm 


4-21 


dependent  on  any  particular  features  of  the  target,  the  adaptive  optics  compensation 
would  have  little  effect  on  the  tilt.  Also,  the  movement  of  the  x  tilt  error  and  y  tilt 
error  as  a  result  of  the  beacon  movement  was  very  similar  to  that  of  the  point  source 
results.  The  x  tilt  error  tended  to  move  away  from  the  adaptive  optics  beacon  in 
the  same  fashion  in  the  point  source  results.  This  could  be  primarily  a  result  of  the 
method  of  tilt  removal  before  the  adaptive  optics  compensation.  The  results  for  the 
y  tilt  error  indicate  the  same  sort  of  bias  as  in  the  point  source  results  as  the  adaptive 
optics  beacon  is  moved.  This  does  not  necessarily  indicate  a  degradation  in  tracking 
performance  (if  it  did  the  tilt  variance  would  have  increased  also)  but  does  indicate 
the  movement  of  the  adaptive  optics  beacon  has  an  effect  on  the  mean  placement 
of  the  scoring  beam  on  the  target  board.  There  is  nothing  in  the  simulation  setup 
itself  that  moves  in  the  y  axis,  therefore  the  bias  must  be  a  result  of  the  movement 
of  the  beacon. 

4-5  Summary 

This  chapter  looked  at  the  results  from  the  simulation.  The  most  notable 
feature  is  the  slope  of  the  tilt  error.  The  error  for  the  X  axis  tends  to  the  positive 
as  the  separation  increases.  The  adaptive  optics  beacon  however  is  actually  moved 
in  the  negative  direction.  The  Y  axis,  on  the  other  hand,  does  not  appear  to  be 
affected  by  the  displacement  of  the  adaptive  optics  beacon  at  all.  This  indicates 
that  the  trend  for  the  X  axis  likely  had  more  to  do  with  the  system  setup  than  the 
anisoplanatic  compensation,  otherwise  the  Y  axis  would  have  been  affected  in  the 
same  way. 
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Figure  4.18.  Tilt  variance.  Each  sub-figure  shows  the  tilt  variance  for  compen¬ 
sated  and  uncompensated  (a  horizontal  line  since  it  is  independent  of 
beacon  separation )  as  a  function  of  separation  between  the  adaptive 
optics  beacon  and  the  track  beacon.  Beacon  separation  is  measured 
in  values  of  90.  (a)  r0  =  10cm  (b)  r0  =  9cm  (c)  r0  =  8cm  (d)  r0  = 
7cm 
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Compensated  X  Tilt  Error  at  Target  rQ  =  0.1 


Compensated  X  Tilt  Error  at  Target  rQ  =  0.09 


Figure  4.20.  X  Tilt  error.  Each  sub-figure  shows  the  tilt  error  for  compensated  and 
uncompensated  (a  horizontal  line  since  it  is  independent  of  beacon 
separation )  as  a  function  of  separation  between  the  adaptive  optics 
beacon  and  the  track  beacon.  Beacon  separation  is  measured  in  values 
of  9q.  (a)  r  o  =  10cm  (b)  r  o  =  9cm  (c)  r  o  =  8cm  (d)  r0  =  7cm 
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Figure  4.22.  Y  Tilt  error.  Each  sub-figure  shows  the  tilt  error  for  compensated  and 
uncompensated  (a  horizontal  line  since  it  is  independent  of  beacon 
separation )  as  a  function  of  separation  between  the  adaptive  optics 
beacon  and  the  track  beacon.  Beacon  separation  is  measured  in  values 
of  d0  ■  (&)  fo  =  10cm  (b)  r  o  =  9cm  (c)  r0  =  8cm  (d)  r0  =  7cm 
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V.  Discussion  and  Future  Work 


5. 1  Discussion 

There  are  a  couple  of  interesting  ideas  to  come  from  this  thesis.  First,  the 
development  of  a  technique  to  measure  the  effect  of  anisoplanatic  compensation  on 
tracking  was  developed.  This  method  allows  for  the  further  analysis  of  other  types  of 
scenarios  to  be  explored.  Second,  and  most  importantly,  it  appears  that  there  is  no 
detrimental  effect  to  the  performance  of  a  centroid  tracker  by  compensating  higher 
order  turbulence  from  an  uncorrelated  angle.  This  is  confirmed  both  by  the  results 
from  the  point  source  data,  and  the  extended  source  data.  This  seems  to  confirm 
previous  work  [13]  that  shows  there  is  very  little  correlation  between  tilt  (Zernikes  Z 2 
and  Z%)  and  the  other  Zernike  polynomials.  This  would  explain  why  the  tilt  variance 
did  not  increase  appreciably  as  the  higher  order  adaptive  optics  beacon  moved  away 
from  the  track  target. 

Another  lesson  that  can  be  learned  from  this  research  involves  the  method  of 
determining  the  level  of  turbulence.  The  turbulence  levels  used  for  this  effort  were 
based  on  r0,  which  was  shown  to  be  an  inaccurate  method  for  horizontal  paths.  The 
strong  turbulence  affected  the  performance  of  the  adaptive  optics  compensation.  Ad¬ 
ditionally,  the  intensity  variations  in  the  images  of  the  both  the  track  beacon  and 
target  board  would  potentially  have  a  detrimental  effect  on  the  centroid  measure¬ 
ments. 

5.2  Potential  for  Future  Research 

Several  new  topics  of  interest  have  been  explored  in  this  research.  These  topics 
present  new  ideas  for  future  theses. 

5.2.1  Compensate  the  HEL.  In  order  to  get  the  best  quality  irradiance 
at  the  target  board,  it  would  be  necessary  to  compensate  the  outgoing  beam  for 
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higher  order  distortions.  This  could  be  accomplished  by  sending  the  HEL  through 
another  optical  phase  delay  map  with  the  same  commands  used  by  the  adaptive 
optics  system.  It  would  then  be  necessary  to  send  the  HEL  back  through  the  same 
path  that  was  travelled  by  the  the  adaptive  optics  beacon.  This  could  be  accom¬ 
plished  by  offsetting  the  track  command  by  the  difference  between  the  track  beacon 
and  the  adaptive  optics  beacon,  i.e.  A x/L.  This  would  still  allow  the  pointing  to  be 
controlled  by  the  TILL  measurements  while  compensating  the  outgoing  beam. 

5.2.2  Feature  Track  Algorithms.  The  next  logical  extension  on  the  project 
would  be  to  develop  a  feature  tracker  and  compare  the  results  from  that  to  those 
from  the  centroid  tracker.  Since  an  edge  tracker  is  much  more  sensitive  to  higher 
order  aberrations,  it  is  likely  that  there  would  be  some  pronounced  effect  on  the 
performance. 

5.2.3  Utilizing  for  Optimal  Tracking.  If  a  clear  performance  degradation 
is  experienced  using  other  tracking  techniques,  a  new  design  would  likely  be  desir¬ 
able.  This  would  likely  consist  of  a  hybrid  of  using  a  compensated  path  and  an 
uncompensated  path.  By  trading  off  from  one  to  the  other,  depending  on  the  turbu¬ 
lence  strength  and  the  observed  value  of  90,  an  adaptive  solution  could  be  obtained. 
Even  more  desirably  the  information  from  both  paths  could  be  combined  optimally 
to  produce  the  “best”  solution.  This  would  require  extensive  examination  of  the 
cross-correlation  properties  between  the  two  paths. 
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